RANDOM WALKS AND MIXED VOLUMES OF 
HYPERSIMPLICES 



Abstract. Below is a method for relating a mixed volume computation for 
polytopcs sharing many facet directions to a symmetric random walk. The 
example of permutahedra and particularly hypersimplices is expanded upon. 



1. Introduction 

Below is a method for relating a mixed volume computation for polytopes sharing 
many facet directions to a symmetric random walk. The example of permutahedra 
and particularly hypersimplices is expanded upon. 

2. Setup 

A family of polytopes is obtained by fixing a finite set H to index the possible 
facets, a subspace X of R ff orthogonal to 1 = (1, . . . , 1) to contain the polytopes 
and a convex subset Y of R H orthogonal to and containing 1 to index the polytopes. 

Write R + for the nonnegative real numbers and F(y) = X n (R+ — y) = {x £ 
X\xh > —yuih G H} for the polytope in X indexed by y G Y. Note that F is 
superaffine: for every a £ [0, 1], y,y £ Y there is aF(y) + (1 — a)F(y) C F(ay + 
(1 — a)y) where the left side involves the Minkowski sum of polytopes. Typically 
these polytopes are not the same. Faces of the polytopes F(y) will be obtained in 
two ways. 

For every A C H the face of F(y) supported by A is F A (y) = Ifl ((R+~ A x 
{0} A ) - y ) = { x g X\x h > -y h Vh G H,x a > -y Q Va G A}. Write Y A = {y £ 
Y\F A {y) + 0} = 7ri({(y,a;) G Y x X\x h > -y h Vh £ H,x a > -y Q Va £ A} for the 
polytope indexing polytopes with faces supported by A, F A — Xr\(R H ~ A x {0} A ) = 
{x £ X\x a = OVa G A} and d A = dim(i^). Note that F A (y) C — y a translate 
of F!a- 

Write S'X* for the unit sphere in the dual space to X. For every g £ SX* the face 
of F(y) with support direction g is F g (y) — F(y) n g~ l (m) where to = g{F{y)) = 
max a . gi ?( J/ )g(x) is the height of the polytope F{y) in the g direction. To relate these 
set S(g, y) C H to be maximal with F s ^ g ^ (y) — F g (y) so all coordinates of j/) 
are perpendicular to the affine span of F g (y). Set C A = {g £ SX*\3y £ Y, S(g, y) — 
A} the support directions for faces with support A so that if g £ C A and x £ F A (y) 
then g(F(y)) = g[x). 

The function of interest will be the mixed volume of {-F(yi)}i£j if the y\ are 
points in Y, which is denoted by Vi(F(y.)). Recall that this means that for {A^} 
positive real parameters, Vol|j| XiF(jji)) — 1 y. X J vj where the sum of 

polytopes in R" is the Minkowski sum and Vi(F(y.)) — jjj, is the fully mixed (or 
equivalently, multiplicity free) coefficient. 
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3. Theorems 

Lemma 3.1. If i E I then 

V I (F(y.))= Y, 9(F(yi))Vi- {i} (F g (y.)). 
gesx* 

This is a sum rather than a more general integral since the volume term is nonzero 
for only a finite set of directions g. See, for instance, the integral over the mixed 
surface function in [1]. 

Lemma 3.2. If g G Ca then the restriction of g{F(-)) to Ya is also the restriction 
of a linear function in (R H )* to Ya- 

Proof: If a G [0, 1] and y,y G Ya and x G Fa(i/) and x G Fa(jj) then g(F(ay + 
(1 — a)y)) = g(ax + (1 - a)x) = ag(F(y)) + (1 - a)g(F(y)). If the affine span 
of Ya does not contain this is sufficient and otherwise it follows similarly that 
g(F(ay)) = ag(F(y)). □ 

Lemma 3.3. // there is J C I with dim(^ eJ F{yj)) < \ J\ then Vi(F(y.)) = 0. 
Proof: This is immediate from the definition of mixed volume. □ 

Lemma 3.4. If {Ui}i£i are subsets ofY and for every A C H either (first case) 
every i has Ui n Ya G {0, Ui} or (second case) there is B D A with \{i\Ui C Yb}| > 
ds then the restriction ofVi(F(-)) to f\ i Ui is multilinear so there are Li G (R H )* 
with Vi(F(y.)) = f]j Li{lJi) if every y { G U l . 

Call such a product subset of Y 1 a zone. 
Proof: Fix g G SX*, y. G l\ { Ui and i E I. Set A = S(g, y t ) so that g G Ca- By 
the hypotheses, either Ui C Ya or there is B D A with \{j ^ i\Uj C Yb}\ > ds- 
This is because yi G J/j n Ya ^ and Yb C so if the first case does not hold 
then U l %Y B - 

Thus for each g, i and y. G Y\j-n Uj the function 

W=0(F(y i ))V / _ {i} (F s (y.)) 

is a linear function of ?/i in [7^. In the first case this is the content of Lemma 3.2. 
In the second case (*) = from Lemma 3.3 with J = {j ^ i\Uj C Y B }. The lemma 
applies since every y € Uj C Y B has F g (y) C F B (y) C F B - y so that X^j e j ^(j/j) 
is contained in a translate of Fg and has dimension at most d B < \ J\- O 
Write 7Tj : Y 1 — > y /_ {*l for the projection. 

Definition 3.5. T/ie vector space of irrelevant signed measures on Y 1 is generated 
by signed measures \i with supp(^) C FJ^. Uj some zone and some i so that for 
every measurable R C Y 1 ^^ there is £((71^)* (/^/i)) = G R H where fn(z) = 1 if 
7Tj(z) G R and fniz) = otherwise. 

Lemma 3.6. If is irrelevant then J y eyI Vi(F(y.))d/j,(y.) = 0. 

Proof: This follows from Lemma 3.4: Assume that n is a generator of irrelevant 
measures with i and {Uj} as in Definition 3.5. 

f V7(F(y.))rf M (y.) = f Vj(F(|/.))dM»-) 
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= / \\Lj{yj)dn{y.) = lira V / Q Lj(yj)dfj,{y.) 

= ^ y2T\ L j( r »j) / L t { yi )d^{y.) 

t R «i a j¥i Jy.eRc.xUi 

= lim Vn^(^)^( E ((^)*(/K Q M))) = 0. 

I ^- a J ... 

Here the limits are over partitions {R a } of Ylj^i Uj an d r " € i? Q so that r aj e Uj. 
□ 

Lemma 3.6 will be used to compute mixed volumes as follows. Note that if one 
probability measure on Y 1 is obtained from another by allowing one coordinate to 
perform a symmetric random walk within some zone then the difference between 
them will be an irrelevant signed measure. If y. G Y 1 then Vj(F{y.)) is the integral 
of Vi{F(-)) over the delta probability measure supported at the configuration y. 
of \I\ points in Y. Assume y. is contained in some zone J\ i Ui. Construct a new 
measure by allowing one point yi to move within Ui via some symmetric process 
until it reaches a boundary point of [/,. This new probability measure will differ 
from the original one by an irrelevant signed measure and hence the integral of the 
volume function will be the same desired value. This process can then be continued 
making use of the fact that the points at which the new measure is supported 
will all be contained in zones L| ■ Vj for which dim(Vi) < dim(£7j), while the other 
regions Vj D Uj may have become larger. In useful cases this will terminate with a 
measure for which the volume integral is easy to compute and the work will have 
been transferred to following the random process to get the final measure. 

4. Example: Permutahedra 

Permutahcdra are the convex hull of a single symmetric group orbit of a point 
in R™ with the defining (trivial plus irreducible) representation of the symmetric 
group S„. These can have anywhere from 1 to n — 1 orbits of facets under £„. 

For the above setup fix a positive integer n and R C [n — 1] = {1, . . . , n — 1}. 
Take H = {T C [n]\\T\ e R}, X = {£ je[n] 6 t \ £ 2 x t = 0} C R H where S t (T) - 1 if 
i E T and otherwise. Take Y C M. H to be the simplex which is the convex hull of 
{|^y^r}re_R where 8 r (T) = 1 if r = \T\ and otherwise. The polytopes F(y) will be 

permutahedra with (generically) \R\ orbits of facets. Hypersimplices are the case 
R = {l,n— 1} and typically have two facet orbits. 

The cases (such as hypersimplices) in which \R\ = 2 and hence Y is an interval 
are particularly nice. 

Fix R — {r, s} with r < s and a linear homeomorphism j : [0, 1] — > Y from the 
interval to the configuration space Y with j(0) — j^S r . 

Lemma 4.1. //Ilie[n] [ u »o> u n] C [0, 1]™ satisfies the following four conditions then 
Ylij([ u io,un]) C Y n is a zone and every maximal zone occurs in this way. 

1) {uie} C {p t } with t e {0} U [r, s] U {n} and p t = ^-i^'J^-ij ■ 

2) If Uio < Pt < un with r < t < s then there is some j with Ujo — Uj\ = pt- 

3) Ifu i0 < p r < un then \{j\u j0 = Uji = p t }\ > r. 

4) Ifu i0 <p s < un then \{j\uj Q = Uji = p s }\ >n-s. 
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Proof Sketch: If A C H = ( [ " ] ) U (™) write A r = A n and A s = A n ( [ " ] ). 
If Ya = j(b><?]) then there is some Y S ( gyV ) = {j(p)}- 

If A = S(g, y) and Y A = {y A } = {j{p A )} then there is K C [n] with A r = {R e 
(W)\R C if} and A s = {i? e D if}. Write A K = A. 

Compute that x K e F AK (j(p\K\)) where 
propriatc constants a and &. 

If r < < s then = F Aff (j(p\ K \)). If r = |if| and K' C K then 

zk' € Fa k If |.K1 = s and K' D if then a; K ' € F Ak 

These points span the faces Fa k (j jsti )) allowing the dimension computations: 
d(A K ) = if r < \K | < s while d(A K ) = r - 1 if r = \K | and = n - s - 1 if 

|Jf| = s. 

This suffices to work out the asserted zones. □ 

The random walk procedure with R = {r, s} is thus a walk of n particles on 
the interval [0, 1] in which the particles pt with r < t < s absorb the first par- 
ticle to reach them and are then invisible, while p r and p s each absorb the first 
several particles to reach them. The number of particles is the same as the num- 
ber of possible interior absorptions so that the random walk procedure will end 
with a measure supported on configurations for which some particle is at an end- 
point (which will make the mixed volume 0) or those which are permutations of 
{PrY x llt=r+i{Pt} x {Ps} n s - Call this particular mixed volume V and compute 
it separately. Now any mixed volume is V times the probability that the system 
ends in a state of this last form. 

Sometimes this probability is fast to compute. 

Bunched at one end: For convenience write qt = Pt if t € [r, s] , q t = p r if t < r 
and q t — p s if t > s. 



Lemma 4.2. If for every t, y t < q t then V [n] (F(y.)) = V]J te[n] ^. Similarly, if 
for every t, y t > q t then V [n] (F(y.)) = V U te[n] 



Proof: Consider the sequence of measures obtained by moving the objects in order 
until the tth particle reaches either qt (which occurs with probability ^-) or 0. □ 

Bunched somewhere: 

If R — {r, s} and there are r < a < b < s with every a < t < b having p t G {j/i} 
and every yi G [p a >?>6] then there is a fast algorithm to compute VLi(_F(?/.)) by 
walking one point at a time. While in the previous case there was only one nonzero 
state to keep track of the probability for at each step, there are now around n. A 
similar process will be polynomial time for any fixed number of bunches. 

This is particularly clean for R = {l,n — 1} (hypersimplices) since the points 
Pi = ^ are evenly spaced in [0, 1]. If a e [0, 1] and k £ M write a + x k for a + k 
mod-1 or the fractional part of the sum. 

Lemma 4.3. For every a\, . . . , a n £ [0, 1] there is 



k—n—l 




Example 1: (three dimensional hypersimplices) 

Set n = 4, R = {1, 3} so H = {{!}, . . . , {4}, {1, 2,3},..., {2, 3, 4}}, X is a copy of 
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t 4 onarnoH t — III —I — I — I <L\ r — ( 1 _2 1 1 _1 _1 3 _1\ 

Jfi hpamieu uyxi ^ 4' 4 ' 4 ' 4 ' 4' 4' 4 ' 4 / ' 2 V 4 : 4)4141 41 4 ' 4' 4/' 

ro - fl I _3 I -I 3 _I _lA T , - (i I 1 _3 3 _1 _1 : n mS w L.jl p 

V47 4> 4'4' 4> 4' 4 1 4/' 4 V4> 4> 4' 4'4' 4' 4' 4 J: 111 u ^ w-uiic 

F is the convex hull of y x = (2, 2, 2, 2, 0, 0, 0, 0) and y 3 = (0, 0, 0, 0, 2, 2, 2, 2) in M 8 
with the linear homeomorphism 

= (2i, 2i, 2i, 2i, 2(1 - i), 2(1 - i), 2(1 - i), 2(1 - i)) of [0, 1] to Y. Consider the 
cases in which every 4ui = 4j _1 (?/i) is an integer with 1 < Aui < 4ui + i < 3 and 
write y. — j'((iti, tt2, U3)). There are ten such y.. Note that F(j(0)) and F(j(l)) 
are points and hence any mixed volume involving either one is zero, F(j{\)) and 
are (dual) tetrahedra and F (j(f )) is an octahedron. A direct computation 
gives V — 32\/2- The Lemma 4.2 gives: 

v m mb hm= mi =y\ 
v [3 mhhm=vui=vi 

V[s]U((h h !))) = V lM = V l ( as P er the normalization) 
and switching i with 4 — i: 

v [3 mbim = ykM = yi 

while the only one not fitting Lemma 4.2 is: 
Lemma 4.3 gives: 

o + v m (j(({, |, i))) + V[ 3] (i((|, |, |))) + V[ 3] (j((i, 1 1))) = k 

o + v [3] (j((h h !))) + v m ml I !))) + = v; 
o + v m u((l i, |))) + o + y [3] (i((i, |, |))) = v, 

+ Vf 3] (i((i, |, |))) + Vfa (?'((!, |, !))) + = V" and 

o + ^ [3] (j((3,i!))) + + = ^ 

Example 2: (three dimensional nonhypersimplices) 

Set n = 4 and r = {1, 2} so that =4 + 6 and the transitions occur at p\ = h 
which absorbs one particle and pi = \ which absorbs two with F(j(p 1 )) again a 
simplex and F(j(p 2 )) a cube. 



5. QUESTIONS 



Guess: If R = {1, n— 1} and every nvi € [0,n] is integral then '"y ! V[w-i] 
is congruent modulo n to ±n" vt. 

This guess is clearly true for the cases in Lemma 4.2. For th this only leaves 
±2 -2-2 congruent to 4 mod 4 to check. 



Question: Is there a polynomial time algorithm for mixed volumes of hyper- 
simplices? 
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